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Statistical fluctuations of the nuclear ground state energies are estimated using shell model calcu¬ 
lations in which particles in the valence shells interact through well defined forces, and are coupled 
to an upper shell governed by random 2-body interactions. Induced ground-state energy fluctua¬ 
tions are found to be one order of magnitude smaller than those previously associated with chaotic 
components, in close agreement with independent perturbative estimates based on the spreading 
widths of excited states. 
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The series of resonances found in the scattering of slow 
neutrons on heavy nuclei, with typical excitation ener¬ 
gies of 6 to 10 MeV, have energy-fluctuations described 
by Wigner or Poisson distribution, depending on their 
quantum numbers These fluctuation properties of 
the quantal spectra Q can be described in great detail 
using Random Matrix Theory |^. The universality of 
the laws of level fluctuations was postulated from the 
conjecture that classical chaotic systems have the same 
Wigner-type fluctuations properties as predicted by the 
Gaussian Orthogonal Ensamble (GOE) while clas¬ 
sical integrable systems have Poisson-type distributions 
[3. These results represent the strongest evidence of the 
presence of quantum chaos in nuclear systems nil. On 
the other hand, it has been convincingly shown that the 
power spectrum (Eourier transform squared) of the fluc¬ 
tuations of a chaotic quantum energy spectra is charac¬ 
terized by 1// noise [l^. 

Extending these ideas to the nuclear ground state, and 
challenging preconceptions, it was recently postulated 
that the chaotic motion of the nucleons inside the nu¬ 
cleus induces chaotic fluctuations in the atomic masses 
0 , of the same size of the deviations between experimen¬ 
tal masses and those computed using the Finite Range 
Droplet Model (FRDM) of the order of 0.7 MeV. 
This coexistence between order and chaos in the nu¬ 
clear ground state was interpreted as an inherent limit 
to the accuracy with which nuclear masses can be calcu¬ 
lated [ll|. 

The error distribution of the mass formulas of Moller et 
al. 0 have a conspicuous long range regularity mui 
that manifests itself as a double peak in the distribu¬ 
tion of mass differences This striking non-Gaussian 
distribution was found to be robust under a variety of 
criteria ^3. Detailed analysis of nuclear mass errors 


have been performed, using both the FRDM the 
shell-model inspired mass calculations of Duflo and Zuker 
(DZ)[I1, the Hartree - Fock- Bogoliubov mass calcula¬ 
tions of Goriely et al. ^3, as well as the equations con¬ 
necting the masses of a particular set of neighbo ring nu¬ 
clei known as the Garvey-Kelson (GK) relations |lft[l7| . 
The presence of strong correlations between mass errors 
in neighboring nuclei, calculated under a mean-held ap¬ 
proach, has been clearly exhibited, as well as the exis¬ 
tence of a well defined chaotic signal in its power spec- 
trum, when their correlations are analyzed as time series 
[iM[19| . a similar power spectrum analysis demonstrates 
that the inclusion of additional physical contributions in 
mass calculations, through many-body interactions or lo¬ 
cal information, removes the chaotic signal in the discrep¬ 
ancies between calculated and measured masses [l^ [T^ . 
It has also been shown that the average mass error is 
smaller that 100 keV. It represents an upper bound on 
the influence of the presence of chaos, which could pre¬ 
clude precise deterministic predictions of nuclear masses. 

Another relevant question is whether nuclear masses 
far from stability can be predicted with a similar ac¬ 
curacy. Large-scale shell model calculations shed some 
light on this important issue [^. Their predictive 
power seems robust against long-distance extrapolations. 
Masses of 67 light nuclei in the fp shell, many of 
them unstable, were calculated with an average error of 
215 keV 1^. Errors come mostly from isospin violation, 
and do not increase far from stability. Large-scale shell 
model calculations for the N = 126 Po-Pu isotones 
find errors in the binding energies smaller than 50 keV 
along the whole chain, with no increase for unstable nu¬ 
clei, and “imply a high predictive power for ground-state 
binding energies beyond the experimentally known nu¬ 
clei.” 



2 


Knowing that the upper bound for the average error 
in both the calculability and the predictability of nuclear 
masses is smaller than 100 keV, it is still relevant to gauge 
the influence on the ground state of the “chaotic layers” 
associated with quantum chaos in nuclei. 

In order to understand the nature of the chaotic er¬ 
rors, in 1^ a systematic study of fluctuations in nuclear 
masses was carried out using the shell model restricted 
to a single shell for protons or neutrons. Realistic Hamil¬ 
tonians combined with a small random component were 
employed to study the fluctuations of the nuclear masses, 
and its dependence on the size of the Hilbert space and 
on the width of the two-body interaction. Ground state 
energy fluctuations were calculated for 39 nuclei in the sd 
and p/shells, using a combination of realistic and random 
2-body interactions. Fluctuation widths were found to 
follow an approximate pattern [^, similar to the 

one found using the periodic orbit theory (sj. Although 
these results exhibit the close relationship between the 
A-dependence of the 2-body interactions, the many body 
effects and the fluctuations in the ground state energy, 
both the origin and the size of these fluctuations were 
introduced in an artificial way. 

The statistical fluctuations of the ground state energy 
induced by the stochastic behavior of levels in the vicin¬ 
ity of neutron threshold were investigated in using 
both perturbation theory and supersymmetry. By em¬ 
ploying the spreading width to estimate the mixing of 
the ground state with the states in the shell above, a rel¬ 
atively flat error (for A>50) close to 100 keV was found. 
A more detailed calculation using second order perturba¬ 
tion theory and random 2-body interactions shows that 
the fluctuations decrease with increasing mass number, 
with a functional form 

CTperi = 19 A-4/3^eV. (1) 

In the present letter we report the analysis of ground 
state energy fluctuations induced by random 2-body in¬ 
teractions, through large scale shell model calculations. 
The valence space includes two major shells for protons 
and two for neutrons. Single-particle energies are fixed, 
as well as both the residual interactions between par¬ 
ticles in the lower shells and the interactions inducing 
mixing between shells. These forces are taken from the 
best available parameterizations in the literature. Two 
body interactions in the upper shell are taken randomly 
from a Gaussian distribution, whose width is determined 
from known interactions. Random matrices in the up¬ 
per shell represent the presence of chaos in neutron res¬ 
onances 1^1^. 

To determine the widths of the distributions of the 
2-body residual interaction a statistical analysis of real¬ 
istic interactions employed in shell model calculations, 
restricted to a single shell, was performed. These inter¬ 
actions are Cohen-Kurath for the p shell |2^, Brown- 
Wildenthal for the sd shell |^, the KBS interaction for 
the pf shell 1^. the BonnG interaction for the ads shell 
, and the KLS interaction for the hf shell )3n(l . As 


most of these interactions are fixed for a given shell, with 
no dependence on the mass number, their average width 
was associated with a representative mass number Amp 
at the middle of a given shell. A plot of the interaction 
widths (Tint against the mass number was obtained in this 
way, which is shown in Fig. ^ A simple fit suggest the 
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FIG. 1: Average interaction widths as function of the middle 
shell Amp . 

functional form 

amt (A) = 12.54 A-^-^^MeV. (2) 

To estimate the ground state energy fluctuations, shell 
model calculations were performed using the code An¬ 
toine 1^. Single particle energies were assumed 
fixed, as well as the interactions between particles in 
the lower shells, where the realistic forces optimized for 
shell model calculations were employed. The coupling 
between the lower and upper shells was built from the 
KLS interaction . The interactions in the upper shells 
were assumed random, with values taken from a Gaus¬ 
sian distribution centered at zero, and a width aint(A). 
The opening of the upper shells significantly extends the 
Hilbert space and would imply a renormalization of the 
realistic interactions used in the lower shell. This effect 
is important for detailed spectroscopic calculations, but 
does not influence the energy fluctuations, so we will not 
give further consideration to this problem here. 

Performing full shell model calculation in a space as¬ 
sociated with many valence particles in four major shells 
(2 for protones and 2 for neutrons) is a formidable task. 
The size of the Hilbert space is so huge that in most cases 
they are impracticable even for the best available codes. 
For this reason we performed the calculations for 2 pro¬ 
tons and 2 neutrons outside closed shells with N = Z = 
2, 8, 20, 40. The 4 representative nuclei studied were 
®Be, ^°Ne, ^"‘Ti and ^"‘Mo. The valence space included, 
respectively, the following combinations: p — sd, sd — pf, 
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pf — gds, and gds — hf shells. In the first two cases 1000 
calculations were performed, 500 for the third, and 200 
for the last case, varying in each calculation the numeri¬ 
cal values of the interactions in the upper shells. As the 
energy fluctuations become smaller for the heavier nuclei, 
the smaller number of trials does not affect the quality 
of the results, shown in Fig. |21 The calculated ground 
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FIG. 2: Ground state energy fluctuations for two protons plus 
two neutrons outside a closed core. 

state energy fluctuations in nuclei with two protons and 
two neutrons outside a closed core, cr 2 p- 2 n(A), induced 
by the use of random interactions in the upper shell, are 
shown as points in Fig. [3 They can be parameterized 
with a smooth curve (full line) as 

CT2p-2n(A) = 2.92 (3) 

Notice that the A-dependence of the fluctuations of the 
2-body interaction strengths behaves like as seen 

in Fig. Q] Allowing these fluctuations to be present only 
in the upper shells, and to affect the ground-state energy 
through many-body effects in the diagonalization, leads 
to a dependence in the case of the 2 protons and 

2 neutrons outside a closed core. 

The final goal of this letter is to determine the size of 
the ground-state energy fluctuations for nuclei all along 
the periodic table, most of whom have many valence par¬ 
ticles. For a majority of known isotopes, performing shell 
model calculations including two major shells for pro¬ 
tons and two for neutrons is impossible. To circumvent 
this difficulty, we rely on the relationship between the 
g.s.-energy fluctuations of systems with 2 protons and 2 
neutrons outside the core and those with Zy valence pro¬ 
tons and Ny valence neutrons obtained with statistical 
calculations performed with random interactions in one 
major shell |^. This shell model study was restricted 
to a single shell for protons or neutrons, using a realistic 
Hamiltonian combined with a random component. We 
studied the ground state energy fluctuations of 39 nu¬ 
clei in the sd y pf shells. It was shown that, for a wide 


range of mixing, the fluctuations grow linearly with the 
size of the random component. In this regime, the ratio 
between the fluctuations in nuclei with Zy valence pro¬ 
tons and Ny valence neutrons and those with 2 valence 
protons and 2 valence neutrons ranges between 1 and 7, 
where the largest ratio corresponds to mid-shell nuclei. 
The results for a 10% random mixing within a major shell 
are described in detailed in Fig. 7 and Table II of Ref. 
1^. They are summarized in Fig. |21 In what follows 
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FIG. 3: Mass widths in 20 nuclei in the sd- and fp-shells, 
relative to the 2p-2n width. 


the conservative assumption of a factor of about seven 
between the mid-shell fluctuations and those found for 
2p-2n calculations is employed. 

By applying all the pervious results we conclude that 
the ground-state energy fluctuations induced by random 
interactions in the upper major shell are of the order 

cr(A) w 20 A-i'^^MeV. (4) 

This estimate of the random ground-state energy fluctu¬ 
ations, shown with full squares in Fig. 01 has some in¬ 
teresting features. They are better appreciated by com¬ 
parison to the different results for the mass errors and 
fluctuations presented in Fig. ^ For a mass A = 100 the 
induced chaotic fluctuations are of the order of 40 keV. 
Its mass dependence is very close to A~^l^. The aver¬ 
age r.m.s. mass error in the FRDM shown with full 
circles in F^. 01 were interpreted as a chaotic mass com¬ 
ponent in Q, whose estimated errors are depicted with 
triangles. They are both an order of magnitude larger 
than the present estimates, and decrease more slowly as 
a function of A. The curve with stars presents the per¬ 
turbative calculation of Molinari-Weidenmiiler , with 
an odd dependence in the mass number as calculated in 
this reference. The fluctuations estimated in Ref. us¬ 
ing second order perturbation theory are so close to the 
present results that they are not displayed, since they 
would superimpose. 
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FIG. 4: Comparison of theoretical mass errors in the FRDM 
and in GK calculations with different estimates of ground- 
state energy fluctuations, color on line. 


It is worth mentioning that local mass errors, deduced 
from the Garvey-Kelson relations, have the same aver¬ 
age size 0, and have a very similar mass dependence 
|34| . shown with open squares, as the present estimate of 
chaotic energy fluctuations. A careful analysis would be 
necessary to disentangle the consequences of this coinci¬ 
dence. 

Large shell model calculations including random inter¬ 
actions in the upper shells, which mimic the presence 


of quantum chaos associated with neutron resonances at 
energies between 6 to 10 MeV, have been performed. Af¬ 
ter gauging the average size of the interaction widths 
in different mass regions, diagonalizations were done for 
2 protons and 2 neutrons in very large Hilbert spaces, 
including two major shells for protons and 2 for neu¬ 
trons. A conservative scaling, based on random matrix 
studies in a single major shell, was employed to deduce 
an upper bound for the energy fluctuations at mid-shell. 
This estimate is consistent with the mass errors found 
in large shell model calculations along the N=126 line, 
and with local mass error estimated using the Garvey- 
Kelson relations, all being smaller than 100 keV. The 
mass dependence was fitted, and found to be close to 
A“4/3 . This agrees in both size and functional form with 
the fluctuations deduced independently from second or¬ 
der perturbation theory. Being an order of magnitude 
smaller than the mass errors in the FRDM, associated 
with a chaotic component in |^, we conclude that this 
component cannot be related to nuclear chaos as defined 
by neutron resonances. The fact that this large chaotic 
component, whose existence was confirmed by a power 
spectrum analysis, disappears when additional 2-body 
forces or local mass information are employed, strongly 
suggest that it is nota physical phenomenon but rather a 
characteristic arising from the mean field approximation. 
Investigations along these lines are underway. 
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